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Abstract 

We define a class of transformations of the dependent and independent 
variables in an ordinary difference sclieme. The transformations leave the 
solution set of the system invariant and reduces to a group of contact 
transformations in the continuous limit. We use a simple example to 
show that the class is not empty and that such "contact transformations 
for discrete systems" genuinely exist. 

1 Introduction 

In a recent article [1^ we proposed a definition of a contact transformation for 
an ordinary difference scheme (OAS') of order K 

Ea{{xk}, {Vk}, k ^ n + M,n + M + I, . . . ,n + N) ^ (1.1) 
a = 1,2 K = N-M + 1, n,M,7VeZ, N > M. 

The Lie algebra of the symmetry group of contact transformations in T was 
realized by vector fields of the form 

X=£,nd.^+(i^ndy„ (1.2) 

where the functions ^„ and satisfy the following conditions: 

• They depend on J G Z+ points {xn+k,yn+k), k — 0,1, ... J — 1 with J > 2 
for at least one vector field in the symmetry algebra; 
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• The Lie algebra should be integrable to a Lie group. This implies that all 
coefficients in the J — th prolongation of the vector fields X should depend 
only on {xn+k,yn+k)i < k < J, and not on any further points; 

• In the continuous limit the algebra and the symmetry group of the OAS' 
(|l.ip should reduce to the Lie algebra and Lie group of the ordinary dif- 
ferential equation (ODE) that is the continuous limit of (ll.ip . 

It was shown in [Ij that contact transformations satisfying theses conditions 
do not exist. The theorems proven in [1 are difference analogues of Backlund's 
famous theorem [7] from which it follows that contact transformations can only 
depend on the first derivatives. Indeed, those that depend on higher derivatives 
are prolongations of point transformations |7J [T^ . 

The negative result presented in [I] leaves open the possibility that a more 
general class of transformations exists, that takes solutions of the OAS' (|l.ip 
into solutions, and reduces to contact transformations in the continuous limit. 
The purpose of this article is to show that this is indeed the case. 

In Section 2 we consider a general n — th order ODE and replace it by a 
system of two lower order equations. We show that some point symmetries of the 
system give rise to contact transformations for the original n-th order equation. 
In particular we show that the point symmetries of the system u' = v,v" = 
give rise to the 10 dimensional invariance group of contact transformations of 
the ODE u'" = 0. 

In Section 3 we apply the same approach to difference systems. We discretize 
system u' = v,v" = in a manner that preserves a 7-dimensional subalgebra 
of its Lie point symmetry algebra. We then eliminate the variable v from the 
system and obtain a third order OAS for the variable u{x) allowing a group of 
transformations that in continuous limit include contact transformations. 

Section 4 is devoted to conclusions. We suggest there a less restrictive defi- 
nition of contact transformations for OAS then the one given in [1]. 

2 Contact symmetries of an ODE as a point 
symmetries of a system of ODEs 

Let us consider an m-th order ODE 

u^""^ ^ F{x,u,u',...,u^"'-'^'>), m>3. (2.1) 

It can always be replaced by a system of lower order equations, for instance by 
putting 

u' = v, -yt"^!) =F(a;,u, (2.2) 

With this choice contact symmetries of (j2.1[) are reduced to point symmetries 
of the system (|2.2p . The Lie algebra of point symmetries of the system (|2.2p is 
realized by vector fields of the form 

X = ^{x, u, v)dx + 4'{x, u, v)du + i>{x, u, v)dv (2.3) 
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such that the prolongation pA"-~^''X of X annihilates the system (|2.2p on its 
solution set. 

Since the solution sets of the single equation (|2.ip and the system (j2.2p 
coincide the vector fields (12.31) will also generate symmetry transformations for 
the ODE (|2.ip . Returning to the original variables in (j2.3p we have 

X^^ix,u, u')d, + cj){x, u,u')du + cj)'-'^ {x, u, u')du' . (2.4) 

Thus, if ^ or in (|2.3p depend on v then (12.41) will correspond to the first 
prolongation of a contact transformation for the ODE (|2.ip with 

^l){x,u,v) ^ (l)'^'^\x,u,u'). (2.5) 

Proceeding in this manner we find all point symmetries of the system (j2.2l) and 
some, if not all contact symmetries of the ODE (12.11) . 

Let us now consider the converse problem. Given a finite dimensional Lie 
algebra of vector fields of the form (|2.3p . what is the most general system of 
two second order equations of the form (12. 2p . invariant under the corresponding 
group of point transformations? For simplicity we concentrate on the case of a 
second order system, i.e. n = 3 in (j2.ip and (|2.2p . 

We choose a basis {Xi^ . . . ,Xs} for the given Lie algebra and write the 
second order prolongations of the basis vectors as 

pr^^^x, = i,d., + <^,;a„ + + + ^fd,, + + ,(2.6) 

i = l,...,s. 

To find the differential invariants we must solve the system of quasilinear PDEs 
pr(2)x,F(x,w,u,u',w',w",u") = 0, i = l,...,s. (2.7) 
The system can be written as a system of equations in matrix form as 

Mf/ = 0, ^{F^,F^,F,,FlF'^,F'i.F':). (2.8) 

"Strong invariants" are obtained if the matrix M has maximal rank tm = 
min{N^ 7). The number of such invariants is Ns = 7 — rjv/ and they are invariant 
on the entire jet space. 

The rank of M can be r < tm on some submanifold and then further invari- 
ants, "weak invariants" can exist. The number of all invariants is = 7 — r. 

Let us now consider a specific example for which all contact symmetries are 
obtained from the point symmetries of a lower order system. The equation (|2.ip 
is specified to 

u'" = (2.9) 

and the system (|2.2p is 

u'^v, v" ^0. (2.10) 
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Following a standard procedure [5] we find that the Lie point symmetry algebra 
of (|2.10p is isomorphic to the de Sitter algebra o(3, 2) realized by vector fields 
with a basis: 



X2 = xdu + 

X3 = a, 

X4 — x^du + 2xd, 

X5 = xdx - vdy 



V 



(2.11) 



Xe = udu + vdy 

X-j = 2vd^ + v^du 

Xg — x'^dx + 2xudu + 2ud-u 

Xg = 2{xv — u)dx + xv'^du + v'^dy 

XiQ = 2x{2u~ xv)da; + (4:u'^ - x'^v'^)du + 2v{2u- xv)dy. 

We shall denote the corresponding local Lie group G. 

By construction (I2.1ip is the Lie algebra of point transformations of the 
system (|2.10[) . If we eliminate v from (j2.10[) we return to the single third order 
ODE (|2.9p . Eliminating v in the same way from the vector field (12. lip we obtain 
the first prolongation of the Lie algebra of contact symmetries of (|2.9p in the 
form (P^ . 

Specifically for the vector fields (|2.11l) . we see that ^ and (j) do not depend 
on V for Xi, . . . ,Xq and Xg so they will generate point symmetries for (|2.ip . 
On the other hand Xi,Xq and Xiq turn into contact symmetries of (|2.ip after 
the substitution (|2.5p . 

We mention that this is an alternative way of calculating the group of contact 
symmetries for (|2.9p to the standard one, used for example in Ref [3]. The result 
is of course same. 

Let us now consider the converse problem for the subalgebra Lq — {Xi, . . . , Xj} 
where X7 generates a contact transformation for the ODE ()2.9p . 

The matrix M of (j2.8p has vm ~ 7 so the group corresponding to Lq has no 
strong invariants. 

However for u' = v, v" = we have r(Af) = 5 and on this manifold the 
system of ODEs (|2.10p is invariant (or at least "weakly invariant"). Moreover, 
applying the prolongations of Xg, Xg and Xiq to the system (|2.2I) we verify that 
this system is invariant under the entire dc Sitter group G ~ 0(3,2). To sum 
up, the only invariants of G 0(3, 2) and Go in this case are weak ones, namely 



All subgroups of 5*0(3,2) were classified into conjugacy classes in Ref [1]. 
The subgroup Go C G corresponding to the Lie algebra Lq — {Xi, . . . ,^7} is 
isomorphic to one of the 7 maximal subgroups of 5*0(3, 2) called the "optical 
group" Opt{2, 1) of 2+1 dimensional Minkowski space. The algebra Lq was 
already known to S. Lie [HI [12] as was its 6 dimensional subalgebra Lq = 



Ii — u — V ~ 0, 



and 




(2.12) 
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{Xi, X2, X3, X4, X5, X7}. The algebras L,Lo,Lo are the only finite dimensional 
Lie algebras of contact transformations of a complex plane (other than point 
transformations) [TTl [T^ . 



3 Contact symmetries of an OAS' as point sym- 
metries of a lower order system? 

Let us now try to produce an analogue of contact transformations for a difference 
system in the same manner as we did for the ODE (|2.9|) in Section 2. 

We again start from the subalgebra Lq = {Xi, . . . Xj} and look for a dif- 
ference system that allows Lq as a Lie point symmetry algebra. We use the 
formalism for symmetries of OAS as outlined in [3 [S] . The idea is to con- 
struct a 3 point difference system 

Eai{Xk,Uk,Vk}k=n,n+l,n+2) = 0, a=l,2,3 (3.1) 

which has Lq as its symmetry algebra and reduces to the system (I2.10[) in the 
continuous limit. 

In order to facilitate the continuous limit we use the variables 

„(1) _„(1) 

(1) _ + (2) _ Pn+l ^,^ 

Un, Pn+l - a;„, i_a;„ ' Pn+2 - — V>-'^) 

^n+2 

_ Vr.+1-Vr. „(2) _ o g«+2 grt+l 

instead of x^+k, Vn+k, fc = 0, 1, 2. 

The relevant prolongations of the vector fields in the variables have the form: 

prX = indx^ + 0n9«„ + i^ndv^ + ^i+i^ (i) -|- V'i+l^ (i) 

/'n + l 1n+l 

+ '/'i+29(2) +V'i'|25„(^) +A(i)9,„^, + A(2)a,.„^,. 

l'n + 2 1n + 2 

Here A''''), (/jC^) and are calculated using the results given in I, namely 

AC^) = ^„+fc - e„+fc-i (3.3) 

A^) ^ I, \Ti{k-l) (k) 1 , aTc 

'^\ik = —k — 7 WfeA -9„+fe=;fc — 7 2^WjA 

"-11+ j ^i=i j=i 

where is the total difference operator 

A^F(a:„ ,«„,?;„ , p~^l^ , q^^l ^ , . . . ) = {i^ (x„+ 1 , u„+ 1 , w„+ 1 , p|,^|2 > 9n+2 > • • • ) 



-i^(a;„, u„, u„,p|,^ji, ...)}. (3.4) 
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More specifically we have 

prXi = du„ (3.5) 
prX2 = x„du„ + + a (1) 

prX4 = xldu^ + 2xndv„ + {2xn + h„+i)d(i) +2d(i) +25(2) 

Pn + l 'in + 1 Pn + 2 1n + 2 

prXj = 2u„aa;„ + w^a„„ + 2/i„+iq,*j^|ia;i„^i 

4.9;, ^ ^"+1 + ^"+2 „(2) N c) 

/,(2) a, ,/,(2) 

1 + 2 

where 



f'n + 2 '(71 + 5 



^(2) „ 9„(1) „(2) u „(2) (2) , l/„(2) n2/. \2 , r,/ (1) n2 

<Pn+2 - ^9n+l9n+2"-n+l ~ "•n+lPn+29n+2 + 2'>^«+2H"»+2'' +^(,'?„+ij 

-U9„ , „(2) . r,, (2) (2) 

+2w«g„+2 + 2g„+i(7„+2"'"+2 - ^^n+2P„+29„+2 

+ 2('?i+2)^^n+2^n+l 

_ /, /„(2) ^2 r,. / (2) x2 (2) 

Vn+2 - --rLn+l[qn+2) " 2/l„+2 (q„+2 ) - '^ln+lln+2- 

To calculate invariants of the subgroup corresponding to p.Sp we impose 

prXaF(a;„,u„,w„,p|,^|i,g^^|i,p^^j2,(?i+2>'i"+i;'i"+2) =0, (3.6) 
a = 1,...,7. 

In matrix form p.6p can be written as 

MUd — 0, Uj — {Fx^^,F.u^,Fy^,F(i) ,F(i) ,F(2) ,F(2) ,Fh^ F^^ ). 

fn + l 'in + 1 Pn + 2 1n + 2 

As in the continuous case, to get a sufficient number of invariants we must 
restrict to an invariant manifold on which the rank of M is r(M) — 5, rather 
than r{M) = 7 as in the generic case. The invariant manifold in this case is 
given by 



h = - V,, - ig^i =0, /2 = q% = 0. (3.7) 
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and on this manifold the invariants are 



/i = 0, /2 - 



h-^, h^{K+,f{p%-<^l) (3.8) 

An invariant difference scheme that reduces to the system (j2.2p in the continuous 
hmit hn+k is 

Pn+l -Vn- = 0, gffja = 0, K+2 = c/l„+i (3.9) 

where c is an arbitrary real constant and in particular c — 1 corresponds to a 
uniform lattice. 

The OAS' (|3.9p is invariant under the group Gq. 

Let us now eliminate the variable w„ from the system (|3.9p . Taking the 
discrete derivative of the first equation in (j3.9p and using the second equation 

(2) 

9n+2 = we obtain 

<l^ul^P% (3.10) 

thus /i /2 = implies /4 = in (|3.8p . Taking the discrete derivative of 13.101 
we obtain 

(2) '^{hn+1 + hn+2 + Ki+z) (3) /Q i i ^ 

'?n+2 = ^77 — T ^ PnU (^-H) 

and hence we have 

/-ON 

= 0' '^"+2 = (3.12) 

as a consequence of (|3.9I) . 

The difference system (|3.9p . together with its difference consequence (|3.10p is 
thus invariant under the Lie group Gq. Moreover, the equations Ii — I2 — li = 
are invariant under the entire group G ~ 0(3, 2) however I3 = c is invariant 
only under Go ^ Opt{2, 1). For the system (|3.9|) G and Go are groups of point 
transformations . 

Now let us consider the third order difference scheme 

pi+3 = 0, hn+2 = cK+i. (3.13) 

We can obtain its symmetry algebra from Lq by eliminating u„, q^l]^i and q'j^]^2 
from all the expressions in (|3.5I) . From p.9|) and p.lOp we have 

= pi+i - ^hn+lP^^l2- (3.14) 

To get the actual vector fields of the symmetry algebra of the OAS' p.l2p (as 
opposed to their prolongations) we need to keep only the coefficients of dx„ and 
du„- From (|3.5p we see that {Xi, . . . ,Xq} remain as point transformations for 
the OAS, however X-^ corresponds to a contact transformation 



Xj 



Pnil - ^/in+lPi+2) + (Pnll ~ ^^n+lPi+2) (3.15) 
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The third prolongation of X-j wih have the form 

Pn + l Pn + 2 + 3 t t 

(3.16) 

The coefRcients c/)'''), A^^) in (IXTC)) were calculated using Maple and they are 
too long to reproduce here. The important fact is that we have 

0W = cj,^''\pli+^l„K+,+,), 0<j<k (3.17) 

XC^) = X^'\plil^^„h^+,+,). (3.18) 

Thus the algebra cannot be integrated to Lie group as shown in general in our 
previous article [1]. 

________ 

However, on the solution set of the OAS' p.l3p we have p„^^ = and 
pA^^X'/ p. 161) simplifies to 

I'n + l t'n + 2 t'n + 3 



This can be integrated to give a one parameter group of contact transformations, 
namely 



- Jl) - PnU ~{3) _ (3) 

1 + Ap„+2 

hl% = h^+kil + Xp^'U), fc = l,2,3. 
In the continuous case we obtain 

X = X + \u.j;, U ^ U + ^\{Ux)'^ , Ux=Ux. (3.21) 

It is easy to see that the p.2ip leaves the ODE p.9p invariant and the OAS' 
p.l3p is invariant under the transformation p.20p 



4 Conclusions 

The main conclusion is that the definition of contact symmetries for difference 
schemes given in [1] was too restrictive. It required that the Lie algebra of 
contact transformations be integrable to a Lie group on the entire jet space. 
Let us propose a less restrictive definition 
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Definition 4.1 The vector fields U.2\) where ^„ and (j)n are functions of 
{xn+j,yn+j,j — 1,...J} form a Lie algebra of contact symmetries of OAS 
if they satisfy the following conditions 

pr^^^XEa\E,=E,=0, a = 1,2; (4.1) 

• At least one of the vector fields in the Lie algebra has J > 2; 

• The Lie algebra should be integrable to a Lie group at least on the invariant 
surface defined by the OAS' (QHJ); 

• In the continuous limit the symmetry algebra and Lie group reduce to a 
Lie algebra and group of contact symmetries of the corresponding ODE. 

The algebra p.Sp constructed in Section 3 is the algebra of contact symme- 
tries of the OAS' p. 121) . This provides an example of the fact that Definition 1 
is not empty: such contact symmetries of OAS do exist! 

Some more specific conclusions concerning the example of ODE y'" = can 
be drawn. 

• The ODE (|2.9p has a 10 dimensional symmetry group of contact trans- 
formations G. The OAS p.l2p is only invariant under a 7-dimensional 
subgroup of contact transformations Go C G. The equation p„_j_3 = 
is actually invariant under the larger group G, but the lattice condition 
hn+2 = chn+i is Only invariant under Gq. 

The situation is similar for the second order equation y" = 0. The ODE 
is invariant under the group SL(3,R) of point transformations. The cor- 
responding OAS is invariant under a 6-dimensional subgroup of SL(3, M) 
isomorphic to the group of general affine transformations of [6] . 

• The Lie algebras of point symmetries of the OAS and the ODE are real- 
ized by identical vector fields (with the correspondence (x, y) o (x„,?/„). 
The contact symmetry is however different: the derivative p = u' is not 
replaced by p^^i but by p^+i — 5^npi+2 (^^e p.lSp as opposed to Xy in 

(EUl)). 

• The contact transformation Xj for the OAS (|3.12p involves "second order 

contact", i.e. ^„ and 0„ in (j3.15p depend on p[l]_i and p^f_^2- ^^'^ 
continuous limit hn+i — >■ this reduces to first order contact (p = Ux 
only). This is in agreement with Backlund's theorem stating that contact 
transformations for an ODE arc of most of order 1 [7] . 

In [S] the authors introduced the concept of "internal" and "external" sym- 
metries. External symmetries are defined on the entire jet space, internal ones 
only on the submanifold of the solutions of the equation. External symme- 
tries refer only to strong invariants, internal symmetries also to weak ones. As 
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stressed in [S^ Backlund's theorem \T. actually only applies to external symme- 
tries. 

In [T] we have shown that the only external symmetries of an OAS' are point 
ones. Here we have shown that the OAS' can allow a class of higher symmetries 
that reduces to contact ones in the continuous limit. In the terminology of Ref 
[5] these are internal symmetries. They do not necessarily depend only on first 
order discrete derivatives. 

It remains to determine whether this class of higher symmetries is actually 
useful, in particular whether it can be used to obtain solutions of ordinary 
difference schemes. To answer this question we are planning to study symmetry 
preserving discretizations of nontrivial ODEs that allow symmetry groups of 
genuine contact transformations [13]. 
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